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Abstract 

Modification of the particles in the course of the source evolution is consid- 
ered. Influence of this effect on multiplicities and correlations of the particles 
is displayed, including an enhancement of the production rates and identical 
particle correlations and also back-to-back particle-antiparticle correlations. 
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1 Introduction 



In this paper we consider effects which arise when particles are initially pro- 
duced inside hadronic matter (particle source) and only afterwards, at the 
very last stage of the process, fly away as free particles. Being produced in 
the hadronic medium they represent a part of the medium, rather quasipar- 
ticles than free particles. So their spectrum E(k) changes in the course of 
the source evolution (say expansion). 

Such a picture leads to specific observable effects. First, production rates 
and identical pion correlations (HBT effect) can be amplified significantly. 
Second, back-to-back particle-antiparticle correlations (PAC) appear which 
are not necessarily small if quasiparticles differ essentially from free particles. 

Let us remind in this connection that momentary transformation of quasi- 
particles into free particles was first considered |IJ in the frame of simple 
oscillator model. Chaotic squeezed states of the particles arise in this model 
giving rise to enhancement of single-particle inclusive production and iden- 
tical particle correlations, but antiparticles (and so PAC) are absent in the 
model. On the other hand, the presence of back-to-back PAC was noted 
some time ago for particle production in the vacuum. However in this case 
PAC are suppressed (no source evolution effect in this case). Unsuppressed 
back-to-back correlations arising due to transformation of neutral quasiparti- 
cles into free particles were found recently in ref . |§ . This kind of correlations 
(including charged pions) was also suggested quite recently [|] as a probe of 
disoriented chiral condensate (DCC) formation. Here we consider pion cre- 
ation inside a finite size source which undergoes time evolution (having finite 
lifetime). Both kinds of effects will be found, irrespective of DCC formation. 

Below we consider multiple production processes at high energies. It is 
suggested that after collision an excited volume is formed which undergoes 
evolution. The last (hadronic) stage of the evolution is under consideration. 
So the specific effects discussed below are expected to be absent in those 
scenarios of particle production which do not treat the hadronization stage 
in detail. In particular one can hardly expect these effects in existing event 
generators though they are natural in hydrodynamical model. 



1 



2 Neutral pion production and evolution 

Neutral pions are represented by real valued field </?(x, t). Its decomposition 
has the form 

= / /o_W2/o F M/2 fl k e + a L e I 1 ) 



(2^)3/2(2^)1/2 

where annihilation and creation operators a^, satisfy canonical commuta- 
tion relations. 

To describe the field emission we introduce an effective classical current 
j{x.,t) (it is the practice in HBT effect study 0, [J). The currents encode 
the space-time region of the particle production including time duration of 
the process. Clearly the currents represent some other particles producing 
pions. The current is considered as a random function and subject of averag- 
ing (chaotic source). Production cross-sections and correlations of identical 
particles can be now expressed through the mass shell Fourier transform of 
the currents: 

j(E k ,k) = J dte^Mt) = J d'xe^-^j^t) (2) 
In particular the field correlators are: 

< 4 x «k 2 >o= 2 ^-^- < f{Ei, k 1 )j(E 2 , k 2 ) > (3) 
and, up to small corrections, 

< a kl a k2 > = (4) 

where brackets include statistical averaging and index in Ihs of the equations 
indicates that the medium evolution is not taken into account. 

If production takes place in a medium one has to modify the free prop- 
agation of the fields, say introduce some external potential. In our problem 
the time evolution of the medium is essential, so we confine ourselves to 
the homogenious sources. Then the problem is reduced to propagation of 
the quasiparticles having variable energy, = E^(t). For example one 
may introduce variable mass m(t) changing in the course of evolution with 
m(oo) = m^. 
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Then our model is given by the Hamiltonian 



H 



n 2 {x) + {\Jv) 2 (x) + m 2 {t)v 2 (x) -j(x)<p(x), vr = (p (5) 



This Hamiltonian is similar to that of quantum oscillator with variable fre- 
quency in the presence of external force (which is the current now). Solution 
of the last problem is known for a long time || . It is helpful to represent the 
solution of the equations of motion in the form of canonical transformation 
of the creation and annihilation operators. The only new feature of the field 
theory model given by Eq.5 (in comparison with nonstationary quantum os- 
cillator) is an appearance of two modes k and — k involved in the canonical 
transformation (as it was the case in the original Bogolubov transformation 
from particles to quasiparticles in superfluid ]7|]). 

The resulting Bogolubov transformation solving the model given by Eq.5 
has the form: 




u v \ ( a k (0) + d k (t) 
v* u*J U-k(°) + d-k( t ), 



(6) 



with 



du(t) = 1 1/2 / «(ti)j k (ti) (7) 

where £(t) is the classical solution for the oscillator with variable frequency 
and with initial conditions 



£(0) = 1, e(0) = i#*(0) (8) 

(we suggest that the currents are absent at t < 0) 

Coefficients of the transformation satisfy an equation 

M 2 -M 2 = i (9) 

and can be expressed through the same function £(t). So the form of ^(t) is 
quite essential. In general it is an oscillating function with variable frequency 
and amplitude. What is more important, time reflected wave appears in £(£). 
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In particular keeping E(t) constant in a small interval near some intermediate 
point t = i one gets 

(2Ey 1/2 £(t) = u*e iM + ve~ iM } E = E(t) (10) 

where u, v are Bogolubov coefficients taken at t — t (with their oscillating 
time dependence canceled). 

The appearance of time- reflected wave leads, throught Eq.7 to two types 
of contributions to observable quantities like < a^a >. The first is of the 
form 

J J dhdhe^-^j^jih) (11) 

and the second is 

J J dhdhe^+^mjit,) (12) 

The contributions of the first type lead to unsuppressed correlations which 
depend on the difference E\ — E 2 being maximal at E\ = E 2 . The second 
type contributions depend on the sum E\ + E 2 (see ref. @) and they are 
always suppressed. Below we neglect the second type contributions keeping 
only "large" ones (of the first type). Note that expectation < aa > also 
contains large contributions. 

To get simple explicit results we simplify the above expressions. Firstly, 
the Bogolubov coefficients are taken to be real valued with 

u = coshr, v = sinhr (13) 

Secondly, the solution £(t) is taken at a middle point i, see Eq.10. We also 
suggest that the field expectations are absent at the initial moment t = 
(no quasiparticles up to t = 0; presumably this is the moment of the phase 
transition to the hadronic phase). Then, using Eqs.6,7 we obtain expectation 
values at t = oo: 

(a t (k 1 )a(k 2 )) = cosh(n + r 2 - fx - f 2 ) \ _ (j*(E u \<n)j{E 2 , k 2 )) 

2^/E 1 E 2 

+ sinhri sinhr 2 F(ki — k 2 ) (14) 

(a(ki)a(k 2 )>, = sinh^ +r 2 -f l - f 2 ) == (f(E 1 , kjjfa, -k 2 )> 

2V^2 
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+^sinh(r 1 + r 2 )F(k 1 + k 2 ) (15) 

to be compared with Eqs.3,4. 

In Eqs.14,15 Ti are full evolution parameters and f\ are the parameters 
evaluated for medium evolution from initial moment t = up to a middle 
point t — i , see Eq.10 (roughly f = r/2). The first terms in rhs of Eqs.14,15 
represent contributions of the pions produced by the currents (that is by some 
other particles) which are now enhanced in comparison with no evolution 
production, Eqs.3,4. The second terms in rhs of these equations contain 
form-factor F(k) which is Fourier transform of the initial volume V, 

F(0) = V7(2tt) 3 (16) 

These terms represent the result of decay of the ground state of the medium. 
In another language, they represent parametric excitation of the field oscil- 
lators existing side by side with forced excitation if the oscillator parameters 
depend on time. 

The evolution effect is determined by the evolution parameter r (the effect 
vanishes if r = 0).The value of r for every momentum k depends in general 
on initial energy E in , final energy Ef and characteristic time duration T of 
the hadronic stage of the process. The parameter r is maximal for small time 
duration 

1 ^ (|0 , ET«1 (17) 



and vanishes for ET 3> 1 (adiabatic process). 

As the final result we obtain expressions for inclusive production of neu- 
tral pions. Single-particle distribution is 

l^k {k) = < flt ( k M k » =cosh(2r-2f)i(|j(^k)| 2 )+sinh 2 r^ (18) 

where the first term is the usual production rate (up to pion energy modifi- 
cation) enhanced by cosh factor arising due to source evolution. The second 
term gives an additional contribution due to ground state decay. 

Two-particle inclusive distribution under consideration is given by 

^ d 3 kd 3 k = ( a i a 2 a i a 2) = (alai)(4 a 2) + (ala 2 )(4 a i) + (aI4)( a i a 2) (19) 
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where the expectations in rhs are given by Eqs. 14,15. The first term in rhs 
of Eq.19 is the product of single-particle distributions, the second term gives 
HBT effect and the third term describes back-to-back particle-antiparticle 
correlatios (PAC,see below) arising due to source evolution. 

The relative correlation function which is measured in experiment is given 

by 

The function Chbt0<-1: k 2 ), describing HBT effect, reaches its maximum 
at ki — k 2 = where 

CW(k,k) = l (21) 

as usual though the slope of its momentum difference dependence is modified 
(it is diminished due to r-dependence, that is due to evolution effect, as one 
can see from Eqs. 14, 18). 

The function CpAc(ki,k 2 ) describes PAC effect and gives an additional 
positive contribution to the correlation, this contribution being maximal at 
kx + k 2 = 0. The width of the PAC peak is expected to be close to that of 
HBT peak (compare Eq.14 and Eq.15). Its height depends crucially on evo- 
lution parameter r, vanishing at r = (see Eq.14) that is in the absence of 
the medium effect. The influence of the PAC effect on the function C(ki, k 2 ) 
in Eq.20 at small momentum differences may be noticeable for soft neutral 
pions having small momenta ki,k 2 and consequently small momentum dif- 
ferences ki — k 2 . So one may expect an increase of the correlation function 
C(k!,k 2 ) at all momenta kj for soft neutral pions having |k| < 1/R where 
R is a characteristic size of the source. In particular from Eqs. 14, 15 and 
Eqs. 18-20 we get an estimate: 



C(0,0) « 2 + 



' 1 L - 0) sinhr + sinh2r 2 



ad 3 kj ^ 011111 ' T 2 (2tt) 



(iS) ( k = 0)coshr + ^sinh 2 r 



(22) 



where we put f — r/2 and neglected pion energy modification in the current 
expectations so that the differential cross-sections in Eq.22 having index are 
those without the evolution effect. Analogous equation is valid for tt + — n~ 
correlations (see below). Numerical estimations of the PAC effect for different 
colliding particles will be given elsewhere. 



6 



3 Charged pion production 



Consideration of charged pions is analogous to that of neutral ones. Now the 
field is complex valued and its decomposition reads 



cf>(x) 



d 3 k 



(27T) 3 / 2 (2£ fe )V2 



a k e 



-iEt+ikx , j f iEt—ik-x. 



(23) 



containing annihilation operators a k for particles (say 7r~) and creation oper- 
ators 6 k for antiparticles. Bogolubov transformation has a slightly different 
form in this case, 



a k (t) = u [d k + a k (0)] + v [d\ + bl k (0) 
b k {t) = u \d k + fo k (0)] + v [cf k + a f _ k (0) 



(24) 



containing function o? k which can be found from <i k by the substitution 
j(x,t)-Kf(x,t). 

As a result single-particle distributions of charged particles have the same 
form as those of neutral ones and two-particle distributions take the form 



{a\a\aia 2 ) = {a\ai){a 2 a 2 ) + (a\a 2 ) (<4 a i) 
for identical (like-sign) pions and 

for particle-antiparticle (7r + 7r~) pair with 



{a\a 2 ) = cosh(r! +r 2 -r 1 - r 2 ) 



^={j*(E 1 ,k 1 )j(E 2 M)) 
2^E t E 2 

+ sinh r\ sinh r 2 F(ki — k 2 ) 



(25) 
(26) 



(27) 



(ai& 2 ) = sinh(ri +r 2 — r 1 — r 2 ) 



2\/e 1 E : 



■.(fiE^m,-^)) 



+-sinh(ri + r 2 )F(k 1 + k 2 ) (28) 



(we suggest again that the fields are absent at the initial moment t — 0). 

So like-sign pions have enhanced HBT correlations with additional ground 
state contribution (maximal at k x = k 2 ) and unlike-sign pions acquire back- 
to-back correlations ( maximal at k x = — k 2 ) due to effect of source evolution. 
The terms of the form < a^a 2 > are suppressed for charged pions and they 
were omitted in the above equations. The resulting HBT and PAC contri- 
butions to the correlation functions of the charged pairs have the same form 
as for neutral pions in the previous section. 

The physical interpretation of the above results is clear from the struc- 
ture of the Bogolubov transformation. Quasiparticles in the source partly 
consist of free pairs (particle having momentum k and antiparticle having 
momentum — k ). The pairs have opposite momenta of their constituents 
not to contribute to the total momentum of the quasiparticle (but influence 
its energy). The same is valid for the ground state of the hadronic medium 
(let us remind that we take into account time variation of the homogeneous 
source, say its selfsimilar expansion). The pairs release when the system de- 
cays into free particles. The HBT effect amplification is also dependent on 
these pairs. 

The origin of doubled correlations between neutral pions is also clear. 
Neutral pions are identical particles and simultaneously they are antiparticles 
to themselves. So they show both types of correlations. 

4 Conclusions 

The evolution of the source in the course of the particle production process 
leads to enhancement of the production rates and to corresponding amplifi- 
cation of HBT effect characteristic for identical particles. Moreover the evo- 
lution leads to appearance of back-to-back particle-antiparticle correlations 
(PAC) which are not suppressed in general. The width of the corresponding 
peak (at k x = — k 2 ) is close to that of the HBT peak (at k x = k 2 ) and its 
height is given mainly by the evolution parameter r. In turn the evolution 
parameter (which also determines the enhancement of single-particle distri- 
butions and HBT effect) depends on the spectrum of the quasiparticles in 
the source and on time duration of the production process, being larger for 
soft particles and for fast processes. This parameter can be determined as 
soon as PAC will be found in experiment. 



8 



Let us note in conclusion that the above effects are of rather general 
nature. Field evolution in the expanding Universe ||, Casimir effect in the 
volume with moving boundaries and particle production in strong electric 
field are described in a similar way using Bogolubov transformation. These 
useful analogies are welcomed. 
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